function eq_ss_new = find_new_ss(eq_SS,param,glob,options)

    C1 = .01;
    C2 = 3;

    tt = 0;
    diff = 10;
    
    tmp.D_out = eq_SS.D;
    
    while tt < options.itermax_DC & diff > 1e-5
        tt = tt+1;
        
        param.C = (C1 + C2)/2;
        Cs(tt) = param.C;
        
        % Step 1: Find D that sets 
        tic()
        eq_tmp = solve_D(param, glob, options,tmp.D_out);
        toc()
        
        % Step 2: 
        
        eq_tmp = find_agg_mit(eq_tmp, param, glob, options);
    
        % Update C
        disp('~~ found D ~~')
        disp('~~ Updating C ~~')
        
        figure(47)
        subplot(3,1,3)
        hold off
        plot(Cs); drawnow;
        
        diff = max(abs(param.C - eq_tmp.C))
        disp(eq_tmp.C)
        
        if eq_tmp.C > param.C 
           C1 = param.C; 
        else
           C2 = param.C;
        end
       
        
       format short
       fprintf('D_out was %.3f \n', eq_tmp.D)

    end
    
    eq_ss_new = eq_tmp;
    eq_ss_new.C = eq_tmp.C;
    
    
    
end